Development in modeling cyclic loading of sands based on kinematic hardening

In this paper, there is presented an elastoplastic constitutive model to predict sandy soils behavior under monotonic and cyclic loadings. This model is based on an existing model (Cambou‐Jafari‐Sidoroff) that takes into account deviatoric and isotropic mechanisms of plasticity. The flow rule used in the deviatoric mechanism is non‐associated and a mixed hardening law controls the evolution of the yield surface. In this research the critical state surface and history surface, which separates the virgin and cyclic states in stress space, are defined. Kinematic hardening modulus and stress–dilatancy law for monotonic and cyclic loadings are effectively modified. With taking hardening modulus as a function of deviatoric and volumetric plastic strain and with defining the history surface and stress reversal, the model has the ability to predict the sandy soils' behavior. All of the model parameters have clear physical meanings and can be determined from usual laboratory tests. In order to validate the model, the results of homogeneous tests on Hostun and Toyoura sands are used. The results of validation show a good capability of the proposed model. Copyright © 2009 John Wiley & Sons, Ltd.     

Historical tsunami in the Makran Subduction Zone off the southern coasts of Iran and Pakistan and results of numerical modeling

Tsunami hazard in the Makran Subduction Zone (MSZ), off the southern coasts of Iran and Pakistan, was studied by numerical modeling of historical tsunami in this region. Although the MSZ triggered the second deadliest tsunami in the Indian Ocean, among those known, the tsunami hazard in this region has yet to be analyzed in detail. This paper reports the results of a risk analysis using five scenario events based on the historic records, and identifies a seismic gap area in western Makran off the southern coast of Iran. This is a possible site for a future large earthquake and tsunami. In addition, we performed numerical modeling to explain some ambiguities in the historical reports. Based on the modeling results, we conclude that either the extreme run-up of 12–15 m assigned for the 1945 Makran tsunami in the historical record was produced by a submarine landslide triggered by the parent earthquake, or that these reports are exaggerated. The other possibility could be the generation of the huge run-up heights by large displacements on splay faults. The results of run-up modeling reveal that a large earthquake and tsunami in the MSZ is capable of producing considerable run-up heights in the far field. Therefore, it is possible that the MSZ was the source of the tsunami encountered by a Portuguese fleet in Dabhul in 1524.     

List of forthcoming papers

List of forthcoming papers     

Exploring the applications of 3D proximity analysis in a 3D digital cadastre

Increasing population in urban areas and limitations of suitable lands for developing houses and urban infrastructure have led to the vertical development in cities. However, these developments are managed by a cadastral system which is mainly two-dimensional and cannot efficiently represent Rights, Restrictions, and Responsibilities (RRRs) in complex scenarios. In fact, a three-dimensional cadastre is required for efficiently registering and representing RRRs. In this paper, a 3D proximity analysis was proposed and implemented to determine RRRs and associated easement rights in non-topology-based data structures. This method can be used to investigate the surrounding spaces of a subject apartment unit or storage in a high-rise. The performance of the developed method was evaluated in a large complex high-rise in Tehran, Iran. The results confirmed that the proposed method can correctly identify the neighbor spaces in complex scenarios.     

Prediction of Land Subsidence Under Cyclic Pumping Based on Laboratory and Numerical Simulations

In this study subsidence due to groundwater withdrawal was investigated. Kerman Province in Iran is struggling with land subsidence problem due to extensive groundwater withdrawal mainly for farming. The rate and type of groundwater withdrawal has very important impact on settlement rate. In this research, effective parameters on land subsidence caused by groundwater withdrawal were determined by laboratory tests. Sampling had done up to depth of 300 m mainly with remolded specimens from Shams-abad, Nouq plain in Kerman province. Similar to the field preconsolidation pressure was applied on specimens in the laboratory. Rate of applied stress on prepared specimens was similar to effect of oscillation of groundwater level. In order to model the actual soil behavior in the laboratory, one-dimensional consolidation device (odometer) was adopted for testing. In these tests, the effect of loading caused by seasonal oscillation of groundwater table is considered by means of cyclic loading in the testing which has great effect on rate of settlements. The results of tests show that when the water table level periodically increases and decreases the amount of settlement decrease, comparing with the case when the groundwater table drop to a constant level. In order to predict the further effects of groundwater level oscillation and actual field condition on land subsidence, a finite element model based on Biots’ three-dimensional consolidation theory was developed. After calibration of finite element model with laboratory tests, this model was used for prediction the effect of groundwater level oscillation on actual field conditions.     

Geochemical Characteristics,Origin, and Evolution of Ore‐Forming Fluids of the Khut Copper Skarn Deposit,West of Yazd in Central Iran

The Khut copper skarn deposit is located at about 50 km northwest of Taft City in Yazd province in the middle part of the Urumieh‐Dokhtar magmatic arc. Intrusion of granitoid of Oligocene–Miocene age into carbonate rocks of the Triassic Nayband Formation led to the formation of marble and a calcic skarn. The marble contains high grade Cu mineralization that occurs mainly as open space filling and replacement. Cu‐rich sulfide samples from the mineralized marble are also anomalous in Au, Zn, and Pb. In contrast, the calcic skarn is only weakly anomalous in Cu and W. The calcic skarn is divided into garnet skarn and garnet–pyroxene skarn zones. Paragenetic relationships and microthermometric data from fluid inclusions in garnet and calcite indicate that the compositional evolution of skarn minerals occurred in three main stages as follows. (i) The early prograde stage, which is characterized by Mg‐rich hedenbergite (Hd_53.7Di_42.3–Hd_86.1Di_9.5) with Al‐bearing andradite (69.8–99.5 mol% andradite). The temperature in the early prograde skarn varies from 400 to 500°C at 500 bar. (ii) The late prograde stage is manifested by almost pure andradite (96.2–98.4 mol% andradite). Based on the fluid inclusion data from garnet, fluid temperature and salinity in this stage is estimated to vary from 267 to 361°C and from 10.1 to 21.1 wt% NaCl equivalent, respectively. Pyrrhotite precipitation started during this stage. (iii) The retrograde stage occurs in an exoskarn, which consists of an assemblage of ferro‐actinolite, quartz, calcite, epidote, chlorite, sphalerite, pyrite, and chalcopyrite that partially replaces earlier mineral assemblages under hydrostatic conditions during fracturing of the early skarn. Fluids in calcite yielded lower temperatures (T < 260°C) and fluid salinity declined to ～8 wt% NaCl equivalent. The last stage mineralization in the deposit is supergene weathering/alteration represented by the formation of iron hydroxide, Cu‐carbonate, clay minerals, and calcite. Sulfur isotope data of chalcopyrite (δ³⁴S of +1.4 to +5.2‰) show an igneous sulfur source. Mineralogy and mineral compositions of the prograde assemblage of the Khut skarn are consistent with deposition under intermediately oxidized and slightly lower fS₂ conditions at shallow crustal levels compared with those of other typical Fe‐bearing Cu–Au skarn systems.     

Mixed precision parallel preconditioner and linear solver for compositional reservoir simulation

Reservoir simulation role in value creation and strategic management decisions cannot be over emphasized. Simulation of complex challenging reservoirs with millions of grid blocks especially in compositional mode is very time-consuming even with fast modern computers. On the other hand, high price of cluster supercomputers prevents them for being commonly used for fast simulation of such reservoirs. In recent years, the development of many-core processors like cell processors, DSPs, and graphical processing units (GPUs) has provided a very cost-effective hardware platform for fast computational operations. However, programming for such processors is much more difficult than conventional CPUs, and new parallel algorithm design and special parallel implementation methods are needed. Using the computational power of CPUs, GPUs, and/or any other processing unit, Open Computing Language (OpenCL) provides a framework for programming for heterogeneous platforms. In this paper, OpenCL is used to employ the computational power of a GPU to build a preconditioner and solve the linear system arising from compositional formulation of multiphase flow in porous media. The proposed parallel preconditioner is proved to be quite effective, even in heterogeneous porous media. Using data-parallel modules on GPU, the preconditioner/solver runtime reduced at least 1 order of magnitude compared to their serial implementation on CPU.     

Optimal water and waste load allocation in reservoir–river systems: a case study

In this paper, a new methodology is developed for optimization of water and waste load allocation in reservoir–river systems considering the existing uncertainties in reservoir inflow, waste loads and water demands. A stochastic dynamic programming (SDP) model is used to optimize reservoir operation considering the inflow uncertainty, and another model called PSO-SA is developed and linked with the SDP model for optimizing water and waste load allocation in downstream river. In the PSO-SA model, a particle swarm optimization technique with a dynamic penalty function for handling the constraints is used to optimize water and waste load allocation policies. Also, a simulated annealing technique is utilized for determining the upper and lower bounds of constraints and objective function considering the existing uncertainties. As the proposed water and waste load allocation model has a considerable run-time, some powerful soft computing techniques, namely, Regression tree Induction (named M5P), fuzzy K-nearest neighbor, Bayesian network, support vector regression and an adaptive neuro-fuzzy inference system, are trained and validated using the results of the proposed methodology to develop real-time water and waste load allocation rules. To examine the efficiency and applicability of the methodology, it is applied to the Dez reservoir–river system in the south-western part of Iran.     

Shigar valley gemstones,their chemical composition and origin,Skardu, Gilgit-Baltistan,Pakistan

A variety of gemstones is being mined in the Shigar valley, Skardu, Pakistan. These include beryl (goshenite and aquamarine), tourmaline (schorl), garnet (almandine–spessartine), apatite, topaz, fluorite, zoisite, clinozoisite, and axinite, mostly occurring in complex or zoned pegmatites and metamorphic rocks. These have been analyzed using electron probe micro-analyzer and X-ray diffractometer. The mineral chemistry of each gemstone is similar to its respective typical gemstone variety with homogenous chemical composition. Field and chemical characteristics suggest that beryl, tourmaline, garnet, apatite, topaz, and fluorite are occurring in zoned pegmatites which are largely formed by magmatic hydrothermal fluids in the cavities and vugs within the intermediate zone. However, zoisite, clinozoisite, and axinite may have a metamorphic and/or metasomatic origin.